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Abstract 

We calculate the low temperature quasi-particle contribution to the ultrasonic attenuation rate in 
the mixed state of d-wave superconductors. Our calculation is performed within the semi-classical 
approximation using quasi-particle energies that are Doppler shifted, with respect to their values in 
the Meissner phase, by the supercurrent associated with the vortices. We find that the attenuation 
at low temperatures and at fields H c \ < H <C H c i has a temperature independent contribution 
which is proportional to y/H where H is the applied magnetic field. We indicate how our result 
in combination with the zero-field result for ultrasonic attenuation can be used to calculate one of 
the parameters vp, H C 2 or £ given the values for any two of them. 
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The discovery of high temperature super conduct ivity9 in the cuprates has led to an intense 
theoretical effort at understanding the origins of the novel phenomena seen in these materials. 
The normal state of the cuprates is highly anomalous! and remains a puzzle that is still 
not understood. However, there is evidence^ that well defined quasi-particles exist in the 
superconducting state. The superconducting state is now known!'! to be a d-wave state 
whose order parameter has d x 2_ y 2 symmetry. The order parameter thus has nodes along the 
lines k x = ±k y in the two-dimensional Brillouin zone. 

While the origin of the pairing interaction that leads to the occurence of d-wave super- 
conductivity is not known at the present time, a great deal of progress in understanding 
the superconducting state can be made by focussing on the low energy nodal quasi-particles 
near the points where the Fermi surface intersects the k x = ±k y lines. These quasi-particles, 
which have an energy dispersion corresponding to Dirac fermionsi in the Meissner phase, 
are the dominant low-energy excitations which determine the low temperature properties in 
the superconducting state. 

It is possible to derive! the d-wave gap from a model Hamiltonian but the connection 
of that Hamiltonian with the underlying physics of strongly interacting fermions is unclear. 
Thus for the purposes of this paper, we assume on phenomenological grounds, the existence 
of BCS like quasi-particles with an energy gap that has d-wave symmetry. These quasi- 
particles might not be the real strongly interacting electrons, to whom their connection is 
unclear at the present time. 

Ultrasonic attenuation has historically been a very useful tool in investigating the energies 
of quasi-particles in conventional superconductors. Early verification! of the BCS prediction 
for the temperature dependence of the energy gap was done using this technique. Even 
afterwards!'0 this has been a useful tool in the investigation of of superconductivity in heavy 
fermion superconductors. On the theoretical front calculations of the ultrasonic attenuation 
in d-wave superconductors have been carried out both in the and dirtyL-aO 

limit in the Meissner phase. The former limit which corresponds to Ql ^> 1 , Q being the 
ultrasound wave-vector and 1 being the electronic mean free path, is the one considered by 
us in this paper. 

Assuming 100 MHz to be a typical frequencyi'0 (although ultrasound experiments can be 
carried out over a large frequency range from kHz to GHz) at which ultrasound experiments 
are done and taking the sound velocity to be§ 4 x 10 5 m/s we find that this translates into the 
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requirement I 3> 4/xm. This restricts the applicability of our work to the cleanest sample 

While ultrasonic attenuation in the vortex state has recieved attention in the past@'@'l§@ 
the emphasis of earlier workers has been to study effects arising from the coupling of ul- 
trasonic waves to vortex motion. This coupling comes from two different effects: a) The 
pinning potentials from impurity ions tend to drag the vortices along with them , b) The 
ionic current due to the ionic displacement will exert a Lorentz force on the vortices. In 
the clean limit (Z > 4/i — m) in which we work, the effects of pinning are expected to be 
weak. The Lorentz force coupling is in any case weak because it is suppressed by a factor 
of v/c where v is the ionic velocity and c is the velocity of light. We therefore ignore these 
processes and focus entirely on the phonon damping from the electronic quasi-particles. 

In conventional s-wave superconductors, low energy electronic excitations in the vortex 
state are the bound states localised in the core region of the vorticesil. However in d-wave 
superconductors, because of the presence of nodes in the gap function, the dominant low 
energy excitations are those in the far region away from the vortex cores. Their contribution 
is expected to be overwhelmingly large than that coming from the the cores. Further, 
experimental studies from STM measurementsBB reveal that that there are just a few 
bound states in the vortex cores. Thus we focus exclusively on the excitations in the far 
region. 

To describe the low energy excitations in the mixed phase we make the use of the semi- 
classical approximation, first discussedi! in the case of s-wave superconductors. This approx- 
imation has been employed for d-wave superconductors in recent timesilS and its regime 
of applicability is the range H c \ < H H C 2, where H is the applied magnetic field and 
H c % and H c2 are the lower critical and upper critical magnetic fields repectively. A strik- 
ing success of this method has been the predictionl! of a term linear in the temperature T 
whose co-efficient scales as y/H (H being the applied magnetic field). This prediction has 
been verified^ and has given greater credibility to the semi-classical description. Recently 
this approximation has been systematized and put on a firm footing by Ramakrishnan and 
Rajagopafi who have derived it microscopically. 

Using the semi-classical approximation to describe the electronic Green functions in the 
superconducting state we calculate the imaginary part of the electron density-density corre- 
lation function which is proportional to the inverse phonon lifetime. As the sound velocity 
is only weakly dependent on temperature^, the temperature dependence of the attenua- 
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tion comes almost entirely from the the inverse lifetime. We find that at low tempeartures, 
and for parameters appropiate to the cuprates, the ultrasonic attenuation co-efficient in 
the mixed state has a temperature independent contribution which scales as vH. The co- 
efficient of this term and its dependence on the ultrasound wave- vector as well as the leading 
temperature corrections to it have been explicitly evaluated. These are our main results. 

We now present the details of our calculation. The semi-classical approximation can 
be understood^ias ignoring, in the first instance, the spatial variation of the supercurrents 
around a vortex outside its core region. This can be justified on the grounds that the spatial 
variation of the electronic wave-function has a characteristic length scale kp 1 which is much 
smaller than the smallest length scale associated with vortex currents £ (the core radius). 
In the cuprates the parameter (fcp£) _1 ~ 10 and thus the semi-classical approximation 
is expected to provide a reasonable description. We then proceed to evaluate the inverse 
phonon lifetime to first order in q where a current has been introduced by taking the energy 
gap function to be A^(Ri, Rj) = A i ^e l ^ Ri+Ri \ Here Ri and Rj are neighbouring sites on a 
square lattice whose fermions are paired up in a singlet state and Ajj = A(— A) for Ri — Rj 
being in the x(y) direction. The spatial dependence of q = (f)/2r is now restored and the 
the inverse phonon lifetime is averaged over a unit-cell of the vortex lattice to get our final 
result for the attenuation co-efficient. 

Solving the Bogulibov-deGennes equation for the mean-field d-wave superconductor to 
linear order in q we obtain the results: 

1 i£ _ v 2 - 

G\{R % - Rj,iuJ n ) = l£ e *(*-*)[ k ~? + |± ] (1) 

where Ujf = (1 + e^/E^)/2, v^ 2 = (1 - e%/E%)/2, eg = £g - /i, £g being the band energy and 



H the chemical potential, E~ = «/e| + A|, Ag = A(cos k x a — cos k y a) being the d-wave gap 
and Ej: = E^ + q- Vfc£g being the Doppler shifted quasi-particle energy. Here 



G^R % - Rj, %uj n ) = -jy r < Q iCT (r)c] i(J (0) > e l ^ T (2) 
is the "normal" Green function. Similarly the anomalous Green functions are given by: 

e iq-(Ri+Rj) e -ik-(Ri~Rj) e ik-{Ri—Rj) 

F ? (R h Rj, iu n ) = Yl u k v k(— p- - — — et) (3) 

and 
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p -iq-(Ri+Rj) _ p -ik-(Ri-Rj) p ik-(Ri-Rj) 

Ff ( R,, R„ i. n) = £ u t , i( — - - —^-) (4) 



k 



Here u^v^ = A^/2E^ and the Green functions F$ and Ft are defined as 

F^Ri, Rj, iu n ) = - [ dr< c id {r)c^(0) > e l ^ T (5) 

J 



F q t(R h Rj, iu n ) = - fdr< 4 T (r) S t a (0) > (6) 

Using the Green functions in equations (1),(3) and (4) we compute the imaginary part 
of the density-density correlation function which is given by 

Xf(Q,u;) = X 1 + X 2 + X 3 + X 4 (7) 

where 

X1 = Z T 5>(^) - n{E U0 )\{ul +e l - ^§^)% + E % - E UQ ) (8) 

k k k+Q 

X2 = = ^ E[»(*U-) - n(E_^)}(vl^vl - ^^)S(u - E_ % + E_^) (9) 

k k k+Q 



and 



^4 = ^ E[l - n(tf s ) - n(F_ fc -_ Q -)](,f + ^| + ^J^W" + E s + E_^). (11) 

k k+Q 

In s-wave superconductors the contribution of X3 and X4 are zero as the ultrasound 
frequency uo <C 2A and so the <5-fn condition in Eqs.(lO) and (11) can never be satisfied. In 
the d-wave case the presence of nodes in Ag means that X 3 and X 4 will be finite. However, 
as the phase space for them is limited, their contribution is expected to be small and so we 
focus exclusively on X 1 and X 2 which make the dominant contribution to x"q(Qi u). 

Now expanding to the leading order in uj we find 
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k k k+Q 

s(e_ s - ^)(4nrf - §§r^ (12) 

k k+Q 

We next expand x to leading order in the ultrasound wave-vector Q and on making the 
substitution k — > — k in the second term in Eq. (12) we arrive at the result 

X?(Q, ") = 2 ^ z n>{El + q ■ V k ^5(Q ■ V k El + QcQp-^-) (13) 

k a p 

where a sum over repeated indices is implied in the arguement of the 5 function. 

To do the k summation in Eq. (13) we only consider! the fermions near the nodes. 
This approximation is reasonable at low temperatures (/cgT <C A). For the node near 
Pa = (ir /2a, 7i /2a) upon linearising the band energies and gap function, we have: £■£ ~ v F k\ 
and A^ ~ v A k 2 where k% and k 2 are co-ordinates normal and tangential respectively to 
the Fermi surface at the node. In terms of the co-ordinates k\ and k 2 we find that the 
contribution from this node to xl/ can De written as 

a2uJ f Jl„ f „/f L.2 7„2 , „.2 7„2 , „ 1 V F k l 



xl a) {QM c^^-jdhj dhnyvlkf + viki + po] 



Vpki -\- v A k 2 



avpkx + (3v A k 2 
o[ — — oi,\Vpk\ — a 2 v F k 2 \ (14) 

lv 2 F kl + v\kl 



Here a = v f Qcos(tx/A — 9), (3 = v A Q sin(7r/4 — 9), a\ = (Qqa 2 /2) cos(ip — 0), «2 = 
— (Qqa 2 /2) cos(ip + 9), p = qv F cos(ip — 7r/4) and 9 and ip are the angles made by Q and q 
respectively with respect to the k\ axis. We now introduce the polar co-ordinates v F ki = 
pcos(j) and v A k 2 = psm<f). Then upon performing the <p integral we find x} reduces to 



x7\Q^) = rdp P n'[p + p }f a [p} (15) 

H TTV F V A J 

where f a [p] = ^ A2 ^ 2)i/2 and A = a-a x p , B = (3 - {y F /v A )a 2 p and p c = ^/nv F v A /a is a 
cutoff introduced to preserve the volume of the Brillouin zone while doing the ^-integration. 
Introducing the variable p' = p + p we have x' a = X + X^ where 



xT\Q^) = -^-rdp'p'ny}f a [p'-p } (16) 



a 2 cu r pc 



' " 71V F V A J po 



and 

xt\Q^) = ~^-Po r dp'n'[p']f a \p' - po] (17) 

1TV F V A J po 

In writing Eqs. (16) and (17) we have made use of p c ^> po which follows from using 
Imax = l/2£ and £ « 30A, -up/a « 192meV and ua/o ~ 28meV^ for parameters appropriate 
to the cupratesi. 

In an analogous fashion we evaluate the corresponding contributions to x'q from the nodes 
near P^ = (— n/2a, — 7t/2a), P c = (—n /2a, it /2a) and Pd = (ir/2a, —ir/2a). We first focus 
on the type of terms in Eq. (16). In this case we find, on expanding to linear order in q, 
that the contribution linear in q exactly vanishes due to cancellations from the terms coming 
from different nodes. We thus find xj" — X$ + + Xg + X<f to be given by 

X « tt^a [ (a 2 + /? 2 ) 3 / 2 (A 2 + r/ 2 ) 3 / 2j 1 ; 

to linear order in q. Once again we have used ksT -C qv F <C p c - Here a and /3 are 

previously defined and we have introduced the parameters A = -upQ sin(7r/4 — 9), r] — 
fA<5cos(7r/4 — 9). The result obtained for this term is identical to the zero field result for 
the attenuation previously calculated in Ref. (12). We next turn our attention to the terms 
of the type written in Eq.(17). Here on adding the contributions from the nodes near P a 
and Pf,, we find to leading order in q, 

>2) "(62) _u^q _ £ gt; F oos(V>-7r/4) 

X,- +X? 7TOA C0S ^ 7F/4j ( a 2 + ^2)3/2 tanll l 2fcflT ) 

Similarly the nodes near P c and Pd yield the contribution 

r > + r ^ co + x/4) ^!^ tanh( £^+iZi) , (20) 

It is straightforward to see that for T — > 0, Eqs. (19) and (20) reduce to 
X^+xT = ~(^) cos^-tt/4) + ^ [g(ooB(^-7T/4))-g(- costy-Tr/4))] (21) 
and 



X^ + X^ = ~(— ) cos(^ + 7r/4) / [g(cos(V»+7r/4)) - g(- cos^ + vr/4))] (22) 

7TC/A ^ i V ) 

Eqs. (21) and (22) are also obtained directly by taking the T = limit for n' in Eq.(15) 
and its analogues for the other 3 nodes. 



To proceed further we now restore q = <p/2r corresponding to the supercurrents around 
a vortex where the vector potential has been ignored as the Ginzburg-Landau parameter 
K ^> 1. Assuming a vortex lattice with circular unit-cells we find that the radius of the cells 



Rc ~ £y H c2 / H. We then average Xq over one un h cell to obtain 



X\H,1)- 7R—T_ r2n ,,,,, ( 2d ) 



to obtain our final result for \ ■ 

In order to perform the averaging in Eq. (23) we make the approximation tanh(x) ~ x 
for \x\ < 1, tanh(x) ~ 1 for x > 1 and tanh(x) ~ — 1 for x < — 1. This approximation 
interpolates between the asymptotically exact behaviour at small and large | x |. Further 
it reproduces the exact result that can be directly calculated for T = 0. In that case the 
integral in Eq. (23) is elementary and can be done by substituting the expressions in Eqs. 
(21) and (22). We then obtain to leading order in the small parameter {^ksT /vf) 2 {H c2 / H)\ 

"(2)/fr T x_ 2uja > 2 i H M/2n k B T 2 H c2 f3 2 

The assumption about the small parameter implies that the regime of validity of our 
result is for kgT ^"("^~) 1//2 together with the condition if cl < H H c2 . For parameters 
relevant to the cuprates H c2 / ' H c \ ~ 100 at low temperatures. This restricts the temperature 
window to T IK at H = H c \. Our result for a temperature independent contribution 
to the attenuation that scales as is entirely understandable as it has its origin in a 
finite density of states at the Fermi-energyH whose size is proportional to \/~H and whose 
signature is seen in the specific heat measurement s@. 

The actual ultrasonic attenuation coefficient is related to \ by 

a s (T,H) = M{Q)x(Q,T,H) (25) 

where M(Q) is a constant which depends on the sound velocity and the electron-phonon 
matrix element and the ultrasound frequency. Then on combining Eqs. (18) and (24) and 
considering the fact that Eq. (18) contains the result for the Meissner phase (H = 0) we 
obtain the result: 

a s (T,H)-a s (T,H = 0) = v F H 1/2 1 _ 8£ 2 ( k B T H c2 

a s (T,H = 0) TiCHj ln2k B T [ 9 1 v F ' H 1 1 ' 



whose window of validity has been described above. This result is remarkable because it is 
independent of the ultrasound wave- vector Q as well as the gap increase parameter v&. Thus 
a measurement of the ultrasonic attenuation at low temperatures and at fields above H& 
would enable a determination of any one of the parameters vp,£ and H c2 given a knowledge 
of the other two. 

We now discuss some shortcomings of our work. We have assumed perfectly well- 
defined quasi-particles thus ignoring the incoherent spectral weight seen in photoemission 
experimentsi. The semi-classical approximation used by us may not accurately describe all 
the physical effects due to the scattering of quasi-particles from the supercurrents in the 
vortex state. Our results are restricted to the clean limit and are applicable only to very 
clean samples. 

Finally we conclude by recapitulating the main points of this paper. We have employed 
the semi-classical approximation to evaluate the phonon damping due to the electronic 
quasi-particles in a d-wave superconductor. We find that for parameters appropiate to 
the cuprates, in a temperature window ksT <C IK, there is a temperature independent 
contribution to the ultrasonic attenuation whose magnitude scales as yH. 
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